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LIST  OF  SYMBOLS 


Ap  Area  of  projectile  base  (m  ) 

C Propellant  charge  volume  (m^) 

E Propellant  gas  covolume  (m'Vkg) 

F Propellant  force  constant  (J/kg) 

G Propellant  gas  ratio  of  specific  heats 

M Projectile  mass  (kg) 

P Chamber  pressure  (MPa) 

PI  Position  of  base  of  projectile  relative  to  breech  face  of  gun  (m) 

PB  Projectile  base  pressure  (MPa) 

3 

R Propellant  density  (kg/m  ) 

SS  Projectile  start  of  motion  (ms) 

Tc  Taylor  coefficient 

U Internal  volume  of  primer  (m-5) 

3 

Vo  Chamber  volume  (m  ) 

Vs  The  number  of  desired  velocity  smoothings  (1,  3 or  9) 

Vc  Cavity  velocity  (m/sec) 

Xp  Projectile  displacement  (m) 

Xc  Cavity  displacement  (m) 

2 

a Projectile  acceleration  (m/sec  ) 

v Projectile  velocity  (m/sec) 

4>  Fraction  of  propellant  burned 
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I.  INTRODUCTION 

In  a continuing  effort  to  reduce  the  cost  and  turn-around  time  for 
the  reduction  of  liquid  propellant  gun  interior  ballistics  data,  a series 
of  computer  programs,  written  in  BASIC  computer  language,  have  been 
developed  to  run  on  a Hewlett-Packard  (HP)  9830A  calculator  with  only 
3808  memory  spaces  available  for  instructions.  The  programs  have  been 
designed  to  reduce  digitized  gun  chamber  pressures  (P)  and  velocities  (v)  ; 
to  compute  projectile  acceleration  (a),  projectile  displacement  (x) , cavity 
penetration  (Xc) , projectile  base  pressure  (Pg) , and  the  fraction-of- 
propellant  burned  (<(>),  and,  to  plot,  scale,  and  label  all  parameters. 

A listing  of  all  programs  is  provided  in  Appendix  A. 


II.  BACKGROUND 

Liquid  propellant  gun  ballistics  data,  in  the  past,  have  been 
reduced  and  plotted  using  the  Ballistic  Research  Laboratories'  Scientific 
Computer  (BRLESC);  however,  several  problems  exist  with  reducing  data 
on  BRLESC  and  among  them  are  cost  and  turn-around  time.  To  address  these 
two  problems,  a series  of  computer  programs  have  been  developed  to  reduce 
and  to  plot  ballistic  data  on  an  HP  9830A  calculator.  Because  of  the 
limited  memory  available  on  the  HP  9830A  calculator,  special  programming 
procedures  had  to  be  developed  to  facilitate  the  desired  data  reduction 
analysis.  The  purpose  of  this  report  is  to  summarize  the  data  reduction 
procedures . 


III.  APPROACH 

Seven  program  files  have  been  created,  on  two  HP  cassette  tapes 
(tape  no.  1 and  tape  no.  2),  to  store  the  programs  required  to  reduce, 
compute,  and  plot  all  desired  parameters.  In  addition  to  the  seven 
program  files  created,  five  data  files  were  created  for  storing  parameter 
which  are  used  for  additional  computations  and/or  for  plotting.  These 
programs  are  designed  to  reduce  liquid  propellant  gun  pressure  and 
velocity  data. 

The  reduction  of  pressure  data  does  not  present  any  particular 
problems,  however,  several  problems  are  encountered  reducing  velocity 
data  and  they  are  as  follows: 

1.  At  firing  pin  contact,  there  is  a voltage  discharge  (across  the 
primer)  pulse  which  is  induced  onto  the  velocity-time  signal  that  may- 
be falsely  read  as  the  start-of-projectile  motion. 

2.  As  the  projectile  approaches  the  muzzle  of  the  gun  there  may  be 
sufficient  gas  leakage  past  the  projectile  that  velocity  measurements 
are  temporarily  lost. 
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These  two  problems,  inherent  in  the  recording  of  projectile  velocity 
data,  are  shown  in  Figure  1.  Figure  1,  a composite  graph  of  a liquid 
propellant  round,  shows  that  at  approximately  0.5  milliseconds  there 
is  an  average  velocity  of  about  300  meters  per  second  and  the 
projectile  has  hardly  begun  to  move.  In  addition,  at  approximately 
3-3.6  milliseconds  (near  the  muzzle),  the  velocity  drops  off  by  a factor 
of  about  two.  It  has  been  assumed  that  neither  of  these  two  events  are 
real.  Therefore,  a scheme  has  been  devised  to  search  for  them  and  to 
linearly  fill  in  the  data  gaps  before  differentiating,  to  compute 
acceleration.  Figure  2 shows  the  results  of  having  replaced  the 
questionable  data  before  differentiating. 

Both  Figures  1 and  2 are  composite  graphs  of  the  same  round  showing 
the  projectile  displacement  (Xp),  and  velocity  (v) ; the  fraction  of 
propellant  burned  (<)>);  and  the  chamber  pressure  (P) . Figure  2 additionally 
shows  the  projectile  acceleration  (a)  and  the  Taylor  cavity  penetration 
(Xc) . The  velocity  data,  shown  in  both  figures,  have  been  smoothed 
nine  times  before  differentiating  to  compute  acceleration.1 

The  fraction  of  propellant  burned  and  the  Taylor  cavity  penetration 
are  based  on  physical  approximation.  The  calculations  of  these  parameters 
are  discussed  in  the  following  section. 


IV.  DATA  REDUCTION  METHODS 

The  numerical  data  reduction  procedures  employed  in  the  various 
programs  are  as  follows: 

A.  Integration 


The  velocity-time  data  are  integrated, 
displacement  (x) , by  using  the  Newton-Cotes 


to  compute  projectile 
quadrature  formula:2 


max 


ydx  = h/3  (vq  + 4v^  + V2) 


(1) 


where  x = 0 at  shot  start  (SS) 
o 

x = total  projectile  travel 
max  r 


^ L . Stansbury , Jr.  and  R.H.  Comer,  "Electronic  and  Numerical  Techniques 
for  Reducing  Noisy  Gun  Ballistic  Data , " BRL  Memorandum  Report  No.  2492, 
June  1979.  (AD  HB006149L) 

2 ■,  ... 

W.  E.  Milne,  "Numerical  Calculus,"  Princeton  University  Press,  Princeton, 

NJ,  P.  123,  (1947). 
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VELOCITY  (M/SEC) 


RACTION 


ACCELERATION  X 1.7  x 10  (M/SEC 

PRESSURE  (MPa) VELOCITY  (M/SEC) 


w 


1 


h = time  increment 
B.  Smoothing 

The  velocity-time  data  are  smoothed  one,  three,  or  nine  times, 
before  differentiation,  using  a moving  arc  polynomial3  over  five 
points : 

73  = (-3 (Vj  ♦ v5)  + 12 (v2  ♦ v4)  ♦ 17(v3))/35  (2) 


where  y^  = the  smoothed  velocity  point,  evaluated  at  the  mid-point  of 
the  five  point  spread. 

C.  Differentiation 

The  velocity-time  data  was  differentiated,  to  compute  projectile 
acceleration,  by  using  the  five  point  numerical  differentiation 

formulas : 4 

y0  12h 


yl'  = 12h 


h4  S 

(-25yQ  + 48y:  - 36y2  + 16y3  - 3y4)  + — y 


(3) 


(-3y0  - 10/j  ♦ 18y2  - 6y3  + y4)  - ^ y5 


(4) 


y 


i 

12h 


(X0  - «yi  + 8y3 


(5) 


where  y'  = first  derivative  of  velocity 

y = number  of  nth  velocity  data  point  (n  = 0 to  4) 

Equations  (3)  and  (4)  are  used  to  compute  the  first  and  second 
acceleration  data  points.  Equation  (5)  is  used  to  compute  all  the 
remaining  acceleration  data  points. 

D.  Interpolation 

The  integration  formula  (1)  is  used  to  compute  the  odd  numbered 
displacement  (x)  points  and  the  following  three  point  Lagrangian  inter- 
polation formula  J is  used  to  determine  the  even  numbered  displacement 
points : 

V.  K.  Sterrett}  Manual  for  Moviri y Polynomial  Aro  Smoothing , BRL  Report 
No.  840 , p.  AS4 , November  1952.  (AD  H109844) 

4 Ibid . 2,  p.  97. 

^T.  B.  Hildebrand,  Introduction  to  Numerical  Analysis,  McGraw-Hill  Book 
Company,  Inc.,  p.  70,  1956. 


x(i)  = °-125x(i-l) 


+ 0 . 75x 


(i+l) 


+ 0.375x 


(i  + 3) 


(6) 


where  i = 2,  4,  6 . . . 

As  the  projectile  approaches  the  muzzle,  there  may  be  a fall-off  in 
the  velocity  signal  due,  in  part,  to  ionized  gases  ahead  of  the  projectile. 
To  bridge  the  data  gap,  the  following  first  order  divided  differences 
interpolating  polynomial6  is  used: 

yi+l  = Xi  ♦ Dl[xi+1  - X]  (7) 


where  i = 0,  1 , . . . m 

and  D1  = ratio  of  first  differences  of  y's  and 
y^  = dependent  variable  preceding  velocity 
x = independent  variable  at  velocity  data 
m = total  number  of  missing  data  points 
E.  Cavity  Penetration 


x' s 

data  fall-off 
fall-off 


The  velocity  (Vc)  and  displacement  (Xc)  of  the  gas  cavity  pene- 
trating into  the  moving  liquid  charge  are  computed  following  the 
procedures  outlined  by  Comer,  Shearer  and  Jones7  and  are  as  follows: 


V 

c 

where  TC 
R1 
a 


TC  yjkl* a 
Taylor  coefficient 
radius  of  apex  of  cavity 
projectile  acceleration 


and,  X = (T*  V At)  + X 
c c p 


where  At 


X 


P 


time  increment 
projectile  displacement 


(8) 


6 Ibid . 2,  p.  210. 

7R.  H.  Comer , R.  N.  Jones  and  R.  B.  Shearer,  "Interior  Ballistic  of 
Liquid  Propellant  Guns,"  BRL  Report  No.  1205,  May  1965  (AD  #344042). 
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F.  Fraction  of  Propellant  Burned 


The  fraction  of  propellant  burned  was  computed  using  the  following 
energy  balance  equation.7 

P/CY-l)  [V  +Ax]  + + C]  v2 

0 pi § (9) 

[(CA/(y-i))  + i^icv  - (PC/Y-1)  (i/p-n) J 

where  P = average  chamber  pressure 

V = internal  volume  of  primer 
o r 

A = area  of  projectile  base 
X = projectile  displacement 
M = mass  of  projectile 
C = mass  of  propellant  charge 

V = projectile  velocity 

Y = ratio  of  specific  heats 

A = propellant  force  constant 
p = propellant  density 
H = propellant  covolume 
$ = fraction  of  propellant  burned 
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5.  Enter  Scratch  A - EXECUTE  into  the  machine  and  rewind  tape. 

6.  Load  file  2,  of  tape  no.  2 and  enter  RUN-EXECUTE  into  machine. 
Program  computes  projectile  displacement  and  stores  data  into  file  4. 

7.  Rewind  tape  and  enter  Scratch  A-EXECUTE  into  the  machine. 

8.  Load  file  3 and  enter  RUN-EXECUTE  into  machine. 

a.  Program  computes  fraction  of  propellant  burned  and  stores  the 
data  into  file  7 of  tape  no.  2. 

b.  See  Appendix  B for  required  input  parameters. 

9.  Enter  Scratch  A-EXECUTE  and  rewind  tape. 

10.  Load  file  6 and  EXECUTE. 

11.  Mount  graph  or  drop-out  blue  paper  on  plotter,  with  pen 
positioned  over  upper-half  of  page,  leaving  enough  margin  for  labeling. 

12.  Enter  Run-EXECUTE.  The  program  smooths  velocity  x-times,  plots 
velocity,  computes  and  plots  acceleration,  and  scales  velocity  and 
acceleration  (V  = 1,  3 or  9) . 

13.  Enter  Scratch  A-EXECUTE. 

14.  Enter  load  file  9 and  EXECUTE. 

15.  Enter  RUN-EXECUTE. 

a.  The  program  computes  and  plots  cavity  penetration. 

b.  See  Appendix  C for  required  input  data. 

16.  Enter  Scratch  A. 

17.  Enter  load  file  5 and  EXECUTE;  however,  do  not  touch  plotter 
until  requested  by  the  machine  (HP) . 

18.  Enter  RUN-EXECUTE.  The  program  (file  5)  does  the  following: 

a.  Plots  and  scales  the  displacement. 

b.  Requests  a "Move  Pen  Down  on  Page"  to  plot  pressure  and 
fraction  burned,  (the  pen  move  is  accomplished  by  pressing  the  "Lower 
left"  plotter  key  board  key,  then  lowering  the  pen  to  approximately 
line  forty-eight  of  the  drop-out  blue) . Ninety  seconds  are  provided 
for  the  move.  This  move  is  required  because  a sheet  of  paper,  on  the 
plotter,  is  divided  into  two  pages  to  accommodate  two  graphs  (x,  v,  a, 

Xc  on  top  graph;  P and  $ on  the  bottom  graph) . 

c.  Plots  and  scales  pressure  and  fraction  burned. 
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d.  Requests  a "Move  Pen  Out  for  P-Labeling."  This  move  is 
accomplished  by  pressing  the  lower  left  plotter  key  and  moving  the  pen 
left  to  the  outer  margin;  then,  press  upper  right  plotter  key  and  move 
pen  to  the  right  margin  (ninety  seconds  are  provided  for  the  move) . 

e.  Requests  a "Move  Pen  Up  for  Labeling."  This  move  is  accomplished 
by  pressing  the  lower  left  key  and  moving  pen  up  to  full  extent  of 
potentiometer. 

NOTE:  For  subsequent  runs,  for  the  same  round,  start  at  item  no.  9. 


VI.  SPECIAL  PROGRAM  NOTES 

As  an  aid  to  the  user,  a series  of  special  program  notes  are  as 
follows : 

1.  The  appropriate  Reel  No.  and  Round  No.  have  to  be  put  into  the 
program  stored  in  file  5,  lines  320  and  330,  respectively,  for  them  to 
be  plotted  onto  the  composite  graphs. 

2.  After  nine  smoothings  and  mid-point  evaluations,  of  the  velocity 
data,  the  first  eighteen  data  points  are  lost;  therefore,  it  is  imperative 
that  the  velocity  reduction  be  started  ahead  of  the  time  of  shot  start 

by  at  least  eighteen  data  points. 

3.  For  the  test  case  (reel  198,  ident  98)  the  velocity  reduction 
was  started  at  the  first  digitized  data  point  of  the  data  block  containing 
the  velocity  data  and  shot-start  was  picked  at  .33  ms  (shot  start  was 
selected  from  the  raw  data  printed  out  as  the  first  data  point  that 
positively  departs  from  the  zero  baseline). 

4.  Time  is  measured  relative  to  the  first  data  point  of  the  second 
block  of  data.  This  was  done  to  prevent  using  additional  memory  space 
searching  for  the  fiducial. 

5.  The  displacement  data  are  multiplied  by  a factor  of  100,  at 
the  time  of  computation,  in  order  that  they  could  be  stored  as  integer 
numbers  and  not  lose  accuracy  at  the  start  of  travel;  however,  the  data 
are  reduced  by  a factor  of  ten  before  plotting. 

6.  Program  changes  are  required  if  new  plotting  scales  are  desired. 

7.  The  acceleration  scale  is  a factor  of  1.7  greater  than  the 
velocity  scale. 

8.  All  plotted  scales  are  labeled  except  the  time  scale. 

9.  The  velocity  data  may  be  smoothed  once,  three  or  nine  times. 

(Each  smoothing  takes  approximately  one  minute  and  fifty  seconds  to  be 
completed) . 
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10.  Only  every  fifth  cavity  penetration  point  is  plotted. 

11.  A search  is  made  of  the  velocity  data  to  determine  at  which 
point  the  data  drops  out  and  when  it  returns.  These  two  prints  are 
saved  and  linear  interpolation  is  used  to  provide  data  points  to  bridge 
the  gap. 


VII.  CONCLUSIONS 

The  series  of  computer  programs,  presented  in  this  report,  have 
been  written  in  the  BASIC  language  and  are  designed  to  reduce  and  to 
display  digitized  liquid  propellant  gun  ballistic  data,  using  an 
HP9830A  calculator  with  only  3808  memory  spaces  available  for  program 
instructions.  These  series  of  programs  are  designed  to  reduce 
ballistic  parameters  (Xp,  v,  a,  P,  <J>,  and  Xc)  previously  reduced  and 
displayed  on  the  BRLESC  computer,  which  required  55K  (72  bits/wordj  memory 
spaces  be  available  for  program  instructions  and  storing  data. 

The  total  time  required  to  reduce  and  display  a composite  graph 
of  desired  parameters,  for  a round,  is  approximately  two  hours  as 
compared  to  forty  minutes  on  the  BRLESC  computer;  however,  the  turn- 
around time,  for  completing  the  reduction  and  display,  on  the  HP9830A 
calculator,  can  be  considerably  less. 

The  reduction  of  liquid  propellant  gun  ballistics  data  on  the 

HP9830A  calculator  offers  two  distinct  advantages  over  reductions  on 
BRLESC  and  they  are,  (1)  reduced  turn-ax-ound  time  and  (2)  reduced 
computational  cost.  However,  due  to  the  limited  memory  spaces,  data 
reduction  on  the  HP9830A  calculator  requires  considerably  more  operator 
interaction,  with  the  calculator,  than  would  otherwise  be  necessary. 
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APPENDIX  A. 


LISTINGS  OF  REDUCTION  AND  PLOTTING  PROGRAMS 


PF5SCEDI  a PaOE  BUNK-NOT  FILMiiD 


REDUCTION  AND  PLOTTING  PROGRAMS 

The  computer  programs,  required  to  reduce  and  plot  liquid  propellant 
gun  ballistic  data,  on  a HP9830A  calculator,  and  the  functions  of  each 
are  as  follows: 

A.  Program  No.  1 (File  4,  Tape  No.  1) 

This  program  reduces  digitized  pressure  (p)  and  velocity  (v)  data 
and  stores  the  data  into  files  0 and  1,  of  tape  no.  2,  respectively. 

B.  Program  No.  2 (File  2,  Tape  No.  2) 

This  program  is  designed  to  search  the  velocity-time  data  for  the 
points  at  which  there  is  data  drop-out  and  when  it  returns.  A second 
search  is  conducted  for  the  points  at  which  the  voltage  discharge  across 
the  primer  interrupts  the  velocity  trace.  The  program  is  further 
designed  to  use  a first  order  divided  difference  interpolating  polynomial 
to  fill  the  data  gaps  before  projectile  displacement  (x)  is  computed. 

The  displacement  data  is  then  stored  into  file  4 of  tape  no.  2. 

C.  Program  No.  3 (File  3,  Tape  No.  2) 

This  program  computes,  "the  fraction-of-propellant  burned"  (4>)  and 
stores  the  data  into  file  7 of  tape  no.  2. 

D.  Program  No.  4 (File  5,  Tape  No.  2) 

This  program  plots  and  scales  displacement  (x) , pressure  (P)  and 
the  fraction  of  propellant  burned  (j>) . The  program  also  labels  x,  a,  v, 

P and  <(>  scales  and  plots  the  round  identification. 

E.  Program  No.  5 (File  o,  Tape  No.  2) 

This  program  smooths  velocity  (1,  3 or  9 times)  and  computes 
acceleration.  In  addition,  the  program  plots  and  scales  the  velocity  and 
acceleration. 

F.  Program  No.  6 (File  9,  Tape  No.  2) 

This  program  computes  and  plots  cavity  penetration. 

G.  Program  No.  7 (File  10,  Tape  No.  2) 

This  program  computes  projectile  base  pressure. 
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PROGRAM  NO.  1 (FILE  4;  TAPE  NO.  1) 


21 


1 f 


y 


1M  DIM  Til  256»2  1 * CC  6 ] * EL  6 D » P I C 256 * 2 ] < V It.  56 « ] 

_ U Q=M= 1 
30  141*0 

40  DISP  THIN-RAW  DflTfl 1 MS ) " * 

50  INPUT  T y 

(.'0  DISP  TMm:  1-h’flW  DflTfl  < MS)  ' ! 

79  INPUT  T9 

U0  DISP  PATH  START  BLOCK1# Vi 
90  INPUT  C9 
100  142=14  3 =G=0 
110  N= 1 
120  2=0 

1 :0  IF  <N  - ?>  THEN  1290 
140  GOTO  160 
150  WAIT  30000 

iso  :;=@ 

170  LOAD  DATA  N»T 

i AO  BC  1 ]=Bt  2 I DC  3 3-BL  4 ]=BI  5 3 ■ B[  6 3=0 
190  J=1 
200  1=1 

Cl]  1+25*  1 H i 1+26*  I3+TE  J+27*  I 3+TI  I 1,  II  [ J 
1 BC  2 3 ;i  1+85*  i I+TC  J+86*  I l+TC  J+87  * I3+TCJ+88*!  +TL  J+  1 , j + | 

1 1 ■ i I I H 161  i MU  J+162*  I3+TI  + 1 i ' +16 

BC  4 3 1 ..  1 l+TC  J+i  26  -I  l+TC  J+  + TC  1 

. 50  1=2 

■ I t 5 > I 3+T  J-t 26*  I1+TC  J+27*  I 

V BC  +•  1-  TC  J+90*  I ]+TC  J+91  * I 3+TC  J+92*  I 3+TC  1+ • , I -TC  l-r--*4-  I 3+TC  +->c-t 
. :0  BC  2 ]=BC  2 ]-BE  1 ] 

. 40  BC  3 3— BC  3 3— BC  1 ] 

00  BC  4 ]=BC 4 3-BC I ] 

10  BC  5 ]=BI  5 3— BC  1 ] 

■20  BC  6 ]=BC  6 3-BC1] 

330  L = 0 
40  B9=BL  l I 
350  Et 1 3=8 

360  IF  < N :•=  5 THEN  40 
370  DISP  MPA  STEP  "* 

::80  GOTO  4O0 

490  DISP  "MP  ■ CAL  STEP"? 

400  INPUT  CC 2 J 
410  FOR  1=1  TO  6 
420  CC I 3=CC 2 3+L 
430  L=L+ 1 
440  NEXT  I 
450  FOR  1=1  TO  6 

460  PRINT  CC I 3 "MPA  OF  MRS" < Bt  I 3"CTS 
470  NEXT  I 
480  PRINT  B9 
490  IF  <C9=2 ''  THEN  510 
300  GOTO  530 
510  N=N+ 1 
520  GOTO  560 
530  IF  (C9=3 > THf  H 550 
540  GOTO  560 
550  N=N+2 
560  LOAD  DATA  II  r 
570  FOR  1=1  TO  3 
580  DISP  ' CAL  L . . i JEFF" I* 

590  INPUT  FC I i 
600  PRINT  FC I I 


1 

i 


k 


22 


■:10  NEXT  l 
■ .0  PRINT 
t30  PRINT 
640  P=Tl=V=B=Iil=U 
-,'ju  IF  < H >=  5)  THEN  700 
’-0  SCALE  0>3j0j8000U  145 

870  hax  i s 0 s l 

1 : 0 .'AXIS  01  10000  145 
890  GOTO  730 
0 O 3 CALE  O * 3 ! O j 8 O 0 0 1 O 
10  :.AXIS  0 j 1 

l AX  I 3 » 1 0 0 0 1 0 , 0 1 0 0 0 ] 1.1 
7 :0  GOTO  740 
40  FOR  1=1  TO  2 
50  FOR  .1=1  TO  256 
60  T 1=T 1+0. 01 
70  TC  Jr  I ]=TC  Jr  I 3 
■-.I  IF  ‘11  . = T 8 • THEN  810 
"90  NEXT  J 
00  NEXT  I 
1 0 PRINT  Jr  I 
830  W2=J 
3:;0  W1  = I 
• 40  i 1 — T 1 
750  FOR  I = U 1 TO  2 
860  FOR  J=W2  TO  256 

IF  CN  - 5 • THEN  40 

0 F [ M ]=  I t 1 ]+  < FI  2 ] + ..  TC  I,  ;;  j 09  >.<  + •;  Ft  3 .♦  It.  -J * I i I ■ 
90  IF  •:  J=25b > THEN  910 
900  GOTO  930 
910  J=1 
920  1-2 
•i  30  GOTO  1000 

4 Vf  0 1 M 1=  < Ft  1 ]+ : Ft  2 3 * 7 T[  J r I 3-B4  1 ■ + • Ft  3 : - Tt  Jr  I ^ -p  - 
0 IF  ' J = 256  ‘ THEN  '470 
" i’  GOTO  1000 


9 7 O J 

1 = 1 

1 =T  1 + 

0 . U 1 

000 

T 1 — T 1 

<*0 . 0 1 

1 0 1 0 

I F •:  T 1 < 0 . 0 1 

7 THEN 

11 

0 

IF  • T1  > = T 

4 1 THEN  1 , 

i 0 30 

IE  'll 

5> 

THEN 

108! 

040 

PLOT 

f i i FT  Q 

j H J 

' * 0 

GOTO 

1070 

; 060 

pr  in  r 

T l j Ft 

0 "MI  - 7' 

j Mr 

1070 

GOTO 

1110 

1080 

PLOT 

T 1 - VC  Q 

j H ] 

1 0 9 0 

GOTO 

1 1 10 

1 100 

PR  INI 

Tl.-VC 

'7 , (1  1 

1110 

IF  •:TI=T9> 

1 HEN  1 

288 

1 120 

IF  • 0 

=256  > 

THEN  1 

140 

1 1 30 

GOTO 

1 160 

1140 

0 = 0 

1 150 

M = 2 

1 160 

GOTO 

1 1 70 

1 1 70 

u=Q+-l 

1 180 

NEXT 

J 

1 1 40 

NEXT 

I 

1200 

PEN 

23 


* 


t U 


J 

i 


I'HEH  1260 


* 


1210  N=N+1 


PROGRAM  NO.  2 (FILE  2;  TAPE  NO.  2) 


10 

DIM  V 

IE  256 <2  1,  > 

at  2! 

20 

C = 0 

3o 

X = 0 

40 

T 1=0 

50 

LORD 

DOTH  1,V 

60 

DISP 

"VEL  DIP'  •' 

j')  5 1 

70 

INPUT 

D9 

80 

IF  CD 

9=1)  THEN 

100 

90 

GOTO 

120 

100  U0SIJE:  990 

110  GOTO  120 

1 20  C=0 . 0000 1*0. 33 3 

1-0  III  P SHOT  3 T hRT  1 MS  ) " 5 

140  INPUT  52 

>50  32=S2  0.011 

• 1 i 1 ES= 1 > 2" 5 

170  INPUT  PI 
o ji  p TMt >:  it  •; 

1J0  INPUT  T9 
20‘  :r  i . i j»0 

-•  I P FPNT-ND  PPOB<?>!  1=Y> 2=N' 
220  INPUT  F9 
. IF  <F9=1  THEN  250 
v 4 * / j i r o 4i  0 

l _ UE;  14';  0 
FOP  1=1  TO  2 
2^0  FOP  .1=1  TO  256 

. I 1=  •'[  ' ’ I 1*10 
_ * IF  > Pl  = l THEN  310 
00  GOTO  320 
'!>  F'P  1 NT  T1 ) VC  J?  I ] 

+0.01 

o If  • i - T9)  THEN  360 
340  NEXT  J 
350  NEXT  I 
360  J=S2 
370  1=1 

1 Dl'-.F'  ' XMFiX  * 1 0 ( M > "5 
>0  INPUT  X9 


400 

T 1 =0+ 

32*0 .01 1 

410 

Q=3 

420 

L=l 

4 30 

XC  Q?  L 

]=XI 0-2, L ]+. l*C VC J, 1 3+4*VC J+l 

440 

IF  <J 

■■=  253)  THEN  50o 

450 

IF  CQ 

=255.'  THEN  560 

460 

IF  (X 

[ 0 j L ] >=  X9)  THEN  629 

470 

J= j+2 

430 

i >=Q+2 

490 

GOTO 

4 30 

500 

i.'  = 0*  j 

510 

XCQ-L 

J=XC  u~2i  l 1+  i.  * VC  25 o > 1 3 + 4* VC  1 

‘3  20 

J=2 

530 

1=2 

540 

Q=Q+2 

550 

GOTO 

430 

56,0 

J=  J+2 

570 

XI  1 > 2 

]=XC  255  ’ 1 1+  \ L * c VC  J*n+4>  a .1+ 1 

530 

J= J + 2 

590 

0=3 

600 

L = 2 

1 '-100 


■ 1 00 


J > - 1 00 


26 


I 


1 O GOTO  430 
620  0=2 
630  L=1 

640  X[  2,  1 3 = -0.  125+XL  5?  1 3+0.75+XC3,  1 3+0.375+XC  1,1  ] 

650  Q=Q+2 

660  XL Q, L ]=-0. 125* X C Q-3, L ]+0. 75+XL 0-1 , L 3+0. 375*XC 0+1 > l 1 
670  IF  (Q+l=255>  THEN  710 
680  IF  <XL  0, L 3 >=  X9-0.25)  THEN  810 
690  0=0+2 
700  GOTO  660 
10  Q=Q+2 

’20  XL 0.  L ]=-0. 125+XL 253,  1 3+0. ?5*XC 255, 1 3+0.375+xC 1,23 
. 30  Q=2 
40  L = 2 

"50  XL  0 , L 3 = - 0 . 1 25 *XC  255 , 1 3+0 . 75*XL  1 , 2 3 + 0 . 8 75+  XL  3 ..  2 j 
760  0=0+2 

0 XL  0,  L 3=-0. 125+ XL  0-3, L 3 + 0. 75  + XC 0-1 , L 3+0. 375-XL Q+l , L 3 
30  U=Q+2 

90  IF  ( XL  0,  L 3 >-  X9-0.25')  THEN  810 
800  G'OTO  770 
8 1 0 s TORE  DATA  4 , X 
820  IF  < P 1 = 1 ■ THEN  840 
30  GOTO  980 
840  IJ=1 
850  Z=S2 

860  FOR  1 = 1 TO  2 . . , • 

870  FOR  J=1  TO  256 
880  PRINT  XL  J, 1 3 ’ T 1 , VL  Z , U 3 
890  IF  < Z = 256 > THEN  '10 
900  GOTO  930 
910  0=2 
929  2=0 
430  T1=T1+0.01 

940  IF  <XC  J, I 3 >=  X9  1 THEN  980 
950  Z=Z+ ! 

960  NEXT  J 
870  NEXT  I 
980  STOP 
990  1=0=1 
1000  .1  = 50 
1010  0=50 
1020  E= 1 

1030  IF  t . A6S 1 VC J+ 1 3, E 3 — V C J ■ I 3 = 87.  7HFN  1 08o 

1040  J=J+1 
1050  0=0+1 

1 060  IF  ( < ,1+18  :>  =2579  THEN  1080 
1070  GOTO  1030 
1030  G= 1 

1090  IF  ( ABS1  VC  G, . 3-  ' ',13*  •=  -'7.  THEN  1180 

1100  G=G+ 1 
1110  0=0+1 
1120  J=J+1 

1130  IF  < J — 2 5 7 • THEN  1150 

140  GOTO  1090 
1150  1=2 
1160  J=1 
1170  GOTO  1090 
1180  L=H=0 

1190  IF  < V£  G,  2 3 VI  -13'  THEN  1,  40 

1200  L=L+ 1 

. 27 


I 


1210 

G=G+ 1 

1220 

IF  03=257)  THEN  1240 

: 2 -'0 

GOTO  1190 

1240 

IF  02  256)  THEN  1260 

1250 

GOTO  1330 

i 260 

1=2 

1270 

Q=Q-256 

1 2 8 0 

H=G-Q 

1290 

V 1 =Y2=B 1 =0 

1 300 

X1=X1=X3=0 

1 3 1 0 

D 1 = B2=N 1 =0 

1 320 

F 1 =F2=0 

1330 

X2=Q-< Q + H 1 

1340 

Y2= <’  VC  0 , I 3-VCG,  I V-/X2  \ 

1350 

PRINT  Vf  r , l ], VC G, I ] 

* 

1 :60 

PRINT  VCt!,  I 3,V2 

0 

1=2 

1 330 

FOR  J=1  TO  H-l 

1390 

V 2 = 22* (. 1 0 + . J <-(0 -- 0 ) ,i 

1400 

B1=V2+V[  (.!,  I ] 

1410 

Vt  Q + J , I :i=Bl 

1420 

PRINT  B 1 , V2 

1430 

NEXT  J 

1440 

.TOT r DAT  ft  l,v 

1450 

RETURN 

1460 

1 = 1 

1 470 

rusp  1 v+5--v=?" ; 

1430 

INPUT  J9 

1490 

K=IJ=S2 

1 500 

IF  ■ VCR +5,  I ]-VC  U j I ] = .J9 ) 

THEN  1 

1510 

K=K+  1 

. 1 >20 

GOTO  1500 

1 5 30 

PRINT  VtK+5,  I 3,VCU,  I ],K 

1540 

M=5+ 1 K-IJ  -+IJ 

1 550 

! 1 = S2-M 

1560 

H1=M-S2 

1570 

:2=M-(M+1 ) 

1580 

F 1 =82  - ( H *■  1 > 

1590 

F2=X 1 -X2 

1 600 

F2,l  1 , X2, XI 

1 6 1 0 

Vis  vr  3.  , I I VC  M,  I 3) /XI 

1 £ 0 

i 2“'  CM,  I i-VC  M + l - I 3>  /X2 

1 620 

n,  vi  n+ 1,13 

i 040 

Iil  = (V  1~V2)/F2 

1 650 

PRINT  IU,Y2-Y1 

1 6 6 0 

FOR  J=1  TO  H 1 - 1 

1670 

+-J>-(S2-0)  ) 

1688 

Bi=Vl  +-VC  82,  I 3 

! f~.  'r1 0 

VCS2+J, I ]=B1 

1 700 

PRINT  B 1 , V 1 

1 7'  1 0 

NEXT  J 

1720 

STORE  Dfllfl  1 i V 

1730 

RETURN 

10  H I M V I L 256  * 2 ] > P I C 256 » 2 3 » X I C 256 » 2 3 » E I C 256 , 2 3 

20  LORD  DATA  0>P 

30  LOAD  DATA  1»V 

40  LOAD  DATA  4>X 

3 0 D I SP  - PR  INTI NG < ? ) - VES= 1 - N0  = 2 " 5 

60  INPUT  PI 

70  D I SP  "PRIMER  VOU  M3) "5 
80  INPUT  U1 


1 00 

INPUT 

V 1 

110 

D I SP 

mog 

120 

INPUT 

M 

130 

D I SP 

’ C < F G ' ' « 

140 

INPUT 

Cl 

150 

D I S F 

GAM  " 1 

16© 

INPUT 

G1 

170 

DISF 

'RHO'PG  M3 

1 

1 80 

INPUT 

R1 

1 99 

DI  P 

AREA  P R 0 J < M 

'i  " 5 

• +J0 

INPUT 

Hi 

210 

DISF’ 

E T A M f G 

" ; 

220 

INPUT 

El 

^ 3 @ 

DISF' 

F < J k G 

34  0 

INPUT 

FI 

250 

s i = i 

PI  • -E 1 

2b0 

D 1 SP 

S MS  1 J 

270 

INPUT 

S*9 

2 2;  0 

89=  <S 

9-0 > 0.01 

2 9 0 

PRINT 

59 

V 0 0 

' w “ 1 . 

. 1 < ■ f:  . 

El  2 = F 2 t 

2 1 0 

PF  IN I 

T i FEE’ V2j F 

iV)  I" 

320 

PRINT 

MS)  ■)  M2S2 ) 

MPA i M S 

3 3 0 

G1=G1 

- 1 

340 

C2=C1 

G1 

350 

C 3 = C 1 

* F 1 Gl 

b 0 

C 4 = ''  1 

R1 : — E 1 

*■  7 0 

1 

Gl 

J 8u 

T 1=0+ 

0 1 0 . 0 1 • +0 . 

01 

*•;  9 0 

2=1 

4 0U 

E=0 

410 

FOR  I 

= 1 TO  2 

420  FOR  J=1  TO  256 

4 ■ 0 :•  2=Pt  30 1 7 b IE + 06  Gl 

■ U ! + ■ R i *Xt  J»  I 3 "100  ■ + C 0. 05i 

:C3+0.0  638*Ci*<<VC S9»?3#ld)t2)“<S2*f  l • 4 
466  EC  J >11=  1 B 1 B.j  1 *■  1 0 0 O 
470  IF  < P 1 = 1 > THEN  490 
480  GOTO  520 

F - 1!  501  1 I f 1 3»  <VI  j . ] ■ 10  ■ , jF  >>Z3 

••  ■ • i 1 ■ 1 F 9 . • 1 1 » F 9 . 0 ) F 9 . t 

510  T 1 =T 1 +0 . u 1 
520  IF  ( 89  = 256  THEN  560 

530  IF  vEC J* I 3 ' *70  THEN  600 

540  39=89+1 
550  GOTO  580 
560  2=2 
570  39=i 
530  NEXT  i 
590  NEXT  I 
600  STORE  DATA  , >t 
610  STOP 

30 


/[  S': 


3*  10  > T-2 


3*101  ■■■' 


PROGRAM  NO.  4 (FILE  5;  TAPE  NO.  2) 


31 


to 
20 
J0 
40 
5 O 
60 


DIM  PIC  256 » 2 ], XII 
LORD  DAT A 0,P 
LORD  DATA  4>X 
LOAD  DATA  7,E 
DISP  "SS(MS) " 5 
INPUT  S8 


256  > 2 ], EIC  £56 


] 


7 O DISP  "TMAXCMS.1 " ! 

80  INPUT  T9 

'<0  DISP  MAX  DISPLACEMENT (M  •"! 


t 00 

INPUT 

X 9 

110 

$9=  $ 

3-4 

3 1 • 0 

.01 

120 

T1=S8 

130 

SCALE 

0 

5 ' 0 

% 4 

140 

XflXiS 

0 

.0.5 

150 

VAX  I S 

er 

1.0 

s 4 

1 6 O 

DISP 

' PRO JE 

tile 

BA 

3E  POS< 

' 0 

INPUT 

PI 

: 0 

FOR  I 

= 1 

TO 

2 

190 

FOR  J 

= 1 

TO 

256 

2 0 0 

PLOT 

ri 

1 * y 

[J , I ] 

l 0 0 ) + P 1 > 

2 1 0 

IF  < X 

j 

1 1 

';=  X9 

■ THEN  250 

220 

T1=T1 

■f  0 

0 1 

S'  so 

NEXT 

j 

240 

NEXT 

i 

250 

PEN 

2 €•  9 

FOR  V 

= 0 

TO 

4 STEP  1 

2 ? 0 

PLOT 

4 . 

4 . Y 

280 

LABEL 

< 

290 . 

1.5,1 

. 7 1 

0 , 10,  3 

290 

FORMAT 

F4.  1 

300 

NEXT 

T 

3 1 0 

DISP 

MOVE 

PEN  DOWN 

ON  PAG 

920 

WAIT 

:'OT0O 

. ::  0 

T1=S8 

340 

SCALE 

0 

. 5 . 9 

, 600 

2 j0 

I PLOT 

O 

J 500 

r,  0 

LABEL 

< 

*■  ) " 

REEL 

198 

" 

7 0 

LABEL 

i’ 

*•  1 

[ D -HO 

50 

:rxis 

( i 

,0.5 

390 

a;  i 

0 

J 100 

400 

FOR  I 

= 1 

TO 

410 

FOR  J 

=s 

9 TC 

256 

420 

PLOT 

T 1 

» PC  . 

, I ] 

4 30 

IF  (T1 

>=  T 

9 1 THEN 

500 

440 

T 1 =T  1 

1-0 

. 01 

450 

IF  CJ 

= £. 

5 6 > 

THEN 

470 

460 

GOTO 

48 

0 

4 70 

39=1 

480 

NEXT 

..1 

490 

NEXT 

1 

500 

PEN 

5 : y 

LABEL  f 

* , l . 

6 1 l . r 

i 0 j 

10/8) 

520 

FOR  V 

-0 

ro 

b 0 0 •: 

— 

m 

— 

1 00 

530 

PLOT 

-1 

j V 

540 

LABEL 

( 

550 

» 

550 

FORMAT 

F 4 , t 

560 

NEXT 

V 

570 

SCALE 

0 

, 5,  e 

, i 

580 

XAXIS 

0 

, 0 . « 

5 9 0 

VAX  I S 

rr 

... 

, 2 / 1 0 

600 

FOR  '-i 

= 0 

TO 

1 STEP  . 

10 

32 


' i PLOT  4.4,  V 
1 11  LRBEL  ( 630  V(' 

1 0 FORMAT  F4.1 

' 40  NEXT  V 
f 00  T1=S8+0.01 
'60  FOR  1-1  TO  2 

■ re  FOR  J=1  TO  256 

' ' PLOT  T1  ,E'[  J,  I 3.-1000 

;,0  IF  < EC  i,  I ] >=  999 > THEN  ’ ;i.i 
’00  Tl=Tl+0.01 
■0  NEXT  l 
. 0 NEXT  I 
'00  REN 

F I'.tR  : • O TO  5 STEP  1 
50  PLOT  X«0.02 
K-0  LRBEL  1 780 >X 
.TO  NEXT  X 
9 FORMAT  F2.0 

■r  '10  I PEN  00  r FOR  f -L  REEL  ] | 

00  WRIT  90000 

10  i RLE  0.5? 0-4 
20  I PLOT  6.5. 1.5 

k;  LRBEl  ■ * < 2 , 1 . 7 j P I 2- 10  8 ! FRACTION  BORNE 

■ K I PLOT  -6 .5)0.  3 

1 • 1 LRBEL  1 * 0. 1 . 7»  p I 2) lO.  8 ' PRESSURE  • Mr  M 
■ Hi  - p MO  ,’t  PEN  UP  FOR  LRBEL  I HO 

:•?£>  WRIT  90000 
• 'TO  SCALE  0 . 5 j 0 . 4 
090  I PLOT  6.5, 1.5 

LRBEL  - • * 1. 7*  PI  2,10  8 • ' DISPLACE  ME- 
1 0 JPH.iT  - 6 . 5 j - 0 . 5 

’•■•0  LRBEL  <*?2j  1 .7»PI/2»  10/8>,'ACCELERATI( 
o.i  I PLOT  O.  S6>  0. 3 

ABEL  • *,  <1.7, PI  2, 10/0' "VELnfTTT  j F 
950  STOP 


in 

DIM  VIC 

256.2  ],  YIC  25 

6 1 2 

2 1.1 

r 1 9 = G 4 = 1 

* 8 

Y=@ 

40 

"IN' 

FOR. 

50 

INPUT  IU 

6 0 

; 10.  MOOTHINGS 

• ( 1 . 

70 

INPUT  38 

80 

DISP  " G 

;£  <.  NS  ;■  5 

,:*0 

INPUT  : 

• l 

100  f 2=1 

1 ‘ 0 . 00  1*B  1 

1 

1 1 

0 DISF' 

PRINTING*?  ■ 

- V E 

12 

0 INPUT 

B9 

13 

0 I SP 

TNi'l.  MS)  ; 

140  INPUT 

Ml 

| cr 

0 lord 

DATA  1 j V 

l 6 

O SCALE 

0,  .11.1800 

'0  MAXIS 

0 .0.5 

:0  YAM  IS 

0.  300 

190  FOR  V 

-o  ro  1800  s 

T EP 

2t 

i0  PLOT 

- 1 . Y 

300 


,'10  LABEL 


. 0 > I 


228 

FORMAT  F 5 . t 

2 3 0 

NEXT  V 

240 

T9=0 . 02 

250 

1 = 1 

268 

.1=02-03=  ;> 

270 

T 1 =0+T9 

288 

YU.  i ]=a 

298 

VC  J.  I 3-YL J 

300 

IF  <T1  >=  • 

GOTO  380 

IF  < G4=S8  > 

338 

GOTO  380 

348 

PLOT  Tl.VC 

350 

IF  ‘IB 9=1  ■ 

368 

GOTO  380 

i | 1(1  1 » 1 

T 1=T 1+U. 01 

590 

J=J+1 

IF  ■ J=255> 

418 

GOTO  280 

420 

M=J-2 

430 

N=  J-l 

440 

U=J+1 

450 

L = J 

460 

1 =0=R=0=B- 

470 

W=2 

488 

YC  I.  I 3=>  l 

490 

VC  J*  I ]=YC 

588 

IF  <G4=S3 

- > ( VC  J-2*  I ]+VC  J+2.  I J * ! •£ 


THEN  340 


■ .<  > VC  M»Q  i VC  1 » W ]>  ‘ 12-  VCN  > P 3 ' Ll » 0 3 > 
13 

THEN  520 

510  GOTO  560 

520  PLOT  Tl.VC Ji 1 3*10 

5 30  IF  CB9=1.'  THEN  550 

540  GOTO  560 

550  PRINT  T1*VC  J - I 3 ■>  J ■■  I 

560  T 1 =T 1 +0 . 0 1 

570  IF  <T1  = W1 ■ THEN  890 

580  M«=M+1 

590  IF  <11=257  THEN  710 
6O0  N=N+1 
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61m 

IF  ■ N=257 

THEN 

740 

620 

0=0+  1 

} ;,(Tl 

IF  » U*25?> 

THEN 

770 

r.  4 0 

l =L  + 1 

*50 

IF  ' L =257 • 

THEN 

800 

t-  6 m 

J=  J + l 

»,?0 

IF  .1=257 

THEN 

8 30 

r£.  Q 0 

1 =1  +1 

r.  *5  0 

IF  •1=257' 

THEN 

860 

’00 

GOTO  480 

710 

M=1 

,'20 

0=2 

30 

GOTO  600 

740 

N=  1 

750 

P=2 

,V.I  i 

GOTO  o>0 

770 

U=1 

780 

0=2 

’9ti 

GOTO  640 

800 

L=  1 

810 

B=2 

: .0 

GOTO  660 

830 

J=1 

840 

I - 2 

850 

GOTO  680 

K=1 

8 70 

W=1 

880 

GOTO  480 

8;  nil 

PEN 

S00 

IF  <C‘9=S8  •• 

THEN 

1 J 00 

9 1 0 

IF  <38=11 

I 1 

1 100 

920 

IF  <38=3 -• 

THEN 

940 

930 

IF  <S8=9 

! HEN 

1010 

940 

02=02+2 

950 

G4=G4+  1 

960 

J=02 

<70 

1 = 1 

T 9=  T 9+0 . 0 2 

' , 9 0 

i 

1000  GOTO  2 7 M 

0 U3=U3+2 

i 02 

0 G4=G4+1 

IO30  j=lj ; 

1040  1=1 

105 

0 T9=T'-'  + m. y 

i 06 

0 DISP  'SMOOTHED 

Q 

10? 

0 WHIT  3000 

! 080  Q9S0--M-1 

! 090  GOTO  2?'i 

1100  DISF  "33 f 

ns  -i 

1110  INPUT  81 

112 

0 IF  <83=11 

THFN 

.0 

! ! ,0  IF  = 

THEN 

80 

140  IF  - 38=3 ■ 

THEN 

0 

115 

0 2=  ■ 3 1-0  • ' 

HI 

116 

0 GOTO  1200 

117 

0 2= <31-01 

HI 

1180  GOTO  1200 

i 190  1-0  ' 

D1 

1200  112  = 2 

— 


l 


1 2 1 0 

RC  1 - 1 3=\ C2*>  - 

25* VC  1 « 1 3+48* VC 

2 i ’ 

if  rt  i - 1 3 o) 

THEN 

1248 

0 

GOTO  1250 

1240 

RC 1- 1 1=6 

1250 

RE  2 * 1 1='  ( £•*'  - 

3*  VC  1 

s 1 3- 10* VC  2 

1760 

IF  C RC  2 j 1 ]< 0 > 

THEN 

1280 

i 70 

GOTO  1290 

1230 

RC  2 - 1 ]=0 

1 2 90 

PRINT  SI +0.O2 

» RC  1 si 

1 3 

1 30O 

PRINT  SI +0.03 

s R C 2 s 

1 3 

1 3 i O 

N-2 

1320 

1 = 1 

1330 

E=1 

1340 

J=1 

1 359 

T 1 =S  1 

1 360 

RC  J+N  i>  I 1-  < C2  + 

<VC  2* 

E 3 -8* VC  2+1 

1 370 

IF  CT1  -=  Ul> 

THEN 

1540 

t-vc  3?  1 3+16*''C  4)13- 


1 1;-!'; 


I : 40  IF  ■ fit J * M ■ I J 0)  THEN  1410 
4O0  GOTO  14£0 
1410  ftt  J +N s I 3=0 


1420 

IF  vZ 

m 

2+1)  THEN 

1430 

IF  <J 

— ^ 

54) 

THEN  1500 

1 440 

Z = Z+ 1 

1450 

J=.J+1 

1 460 

GOTO 

13 

60 

1 470 

E=2 

1480 

2=1 

* 490 

GOT  0 

13 

60 

1500 

J=1 

: 5 1 0 

1=2 

i 520 

N=0 

i 5 30 

GOTO 

1 3 

60 

i‘40 

STORE 

[iftTfl 

8s  R 

0 

T 1=3 1 

+ D 

' Sr.H 

vr  Rt  F 

0 

* 5 » 0 

:0Qm 

i 570 

:r:  ; i 

0 

• 0 . r 

1 530 

VAX  I 

0 

500 

1590 

FOR  1 

= i 

TO 

‘ T-.  Cl  0 

FOR 

= i 

TO 

256 

6 10 

IF  ■£ 

9= 

1 ■ then  1630 

1 620 

GOTO 

it. 

40 

1 630 

PRINT 

T 1 « H[ 

MI-  VC  J . 

1 640 

IF  f Tl 

= S 

1 ■ THEN  1 

1 650 

GOT  0 

It 

70 

1660 

PLOT 

T 1 

> RC  .. 

s I 3 

1.670  IF  • T 1 ■■■=  W1  1 THEN  1710 

1630  T 1 = T 1 + D 1 

1690  NEXT  J 

! 700  NEXT  I 

1710  PEN 

1720  STOP 
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58 


Hi  v..,  ] . ; I [ 256 » 2 3*  CT  250  3 , BC  250  3 

ftp  PftTft  4.  ; 

L HP  DATA  S.H 

0 DI  P TAYLOR  COEF.  5 
50  INPUT  T2 

l-P  C5=0. 001*0. 01 

DIF  ' UftV I TV  RAD I US (M  > "! 

U INPUT  PI 

Dl  F i ,'F £ flPFft''  CM*  • 

00  INF UT  HI 

10  DISF'  "PROJ  BftSE  F'OS.c  M v; 

120  INPUT  B 1 

1 • DIP  SHOT  START • MS'  5 
140  INPUT  Si 

DISP  CAVITY  INITIATION  CM  • '5 
•<■0  INPUT  C8 

-0  D I ,F  PR  INTI NG < ? ■ . 1 =YES > 2 -NO “ ; 

1,0  INPUT  PI 
• ' ?=  i Si  -0  1 0.  01 

200  ft=0 
210  X=0 
220  C9=0 
230  J*S9 
. 4*.'  l=H=i 
50  CC 1 3=BC 1 3=0 
260  0=2 
_ 70  T1=S 1+0.0! 

[ • 5 + + 3QR<  f i *R[  J«  I 3+1000) )+Ct  Q 

0 BI  0 ]•--  CC  0 3+i.  8 + <X£  Mj  I 1 1 
00  r F PI -1  THEN  320 
,10  GOTO  330 

. 0 PIN  I T 1 < C t 0 3-  BC  9 J.  1 ; :t  Qj  I ] +B 1 
30  :[  Q ' I 3=XL  0.  I 3 
40  T 1 =T 1 +0 . 0 1 

50  H=H+1 

IP  • E 3 = '[  0.  13  '108  ■ +F  i THEN  440 

70  J=J+1 
0 4=0+1 

■ -25  1 HEN  410 

40U  GOTO  280 
4lo  1-2 
420  J=1 
430  GOTO  230 
440  . C FiLE  0.5,0.  l 

450  T 1 =S 1 

■ 0 H STEI  5 

4 70  PLOT  Tl.BC I 3 
0 LPLOT  -0.  :.-0„3 
490  LABEL 
500  T l=Tl+0. 05 
510  NEXT  I 
520  PEN 
530  END 
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PROGRAM  NO.  7 (FILE  10;  TAPE  NO.  2) 


f 


L 


40 


A 


lw 
20 
30 
40 
C.  0 
60 
76 
80 
00 

1 00 
1 10 
120 
140 
j 50 
160 
170 
1 80 

1 99 

200 
2 1 0 
220 
;:30 
. '40 

2 50 
260 
270 
280 
290 

; 0 0 
310 
320 

3 30 
340 
356 
360 

: ? O 
380 
390 
400 
410 
420 

4 30 
440 


flU  256.  MIL,  56. 2 3.  PIC  56.2  3 
LORD  DRTfl  0 . P 
LORD  DRTR  7.F 
LORD  DRTR  8.R 
DISP  "PRO, i WEIGHT' G )"! 

INPUT  M 
M=M"453.6 

DISP  "SHOT  START < MS > " < 

INPUT  SI 

DISP  "CHARGE  WT ‘.G  • "! 

INPUT  C 

C = C 4 5 3 . 6 
T 1 =S 1 +0 . O l 
Q=(  SI  -'0. 01 

DISP  "BORE  hPER1  M2)"! 

INPUT  R1 

A1=A1  6.4516E-04 

PRINT  " Ti  BF 1 M > F 

PRINT  MS  MF'R  MFh 

K=0 
E-l 

FOP  1-1  TO  2 
FOR  J=1  TO  256 
BC J . I ]=M*RC  0 . E 3*32 
IF  <1=1. < THEN  320 
Ul  = « M+<  C+.  1 -F[  J.  I ],•  1000  ) 

IF  < FC  .I j I ] 1000  :■•=  0.97) 

GOTO  320 
K=  1 
U 1 -0 

IF  <11  = 4;  THEN  430 

R INI  T 1 ' ET  ,l < I ]*  : 30  145!  Pl  1 1 . E ]j  UM 100  1 4 
Q=Q+ 1 

IF  '-0=257  THEN  370 

GOTO  400 

Q=1 


< 32 . 6 + R 1 > 


BP 


1 + AC  0 ! E i * : 2 
THEN  30O 


5 2 MR  : 


GOTO  230 

Tl=Tl+0. 
ME ' : T J 
NEXT  I 
STOP 
END 


INPUT  PARAMETERS 


1.  Ul,  internal  volume  of  primer,  (m  ) 

2.  VI,  chamber  volume,  (m  ) 

3.  M,  mass  of  the  projectile,  (kg) 

4.  C,  mass  of  the  propellant  charge,  (kg) 

5.  GAM,  propellant  ratio  of  specific  heats 
b.  RHO,  propellant  density,  (kg/m3) 

7.  AREA  PROJ , area  of  base  of  projectile,  (m‘ 

3 

8.  ETA,  propellant  co-volume,  (m  /kg) 

9.  F,  propellant  force  constant,  (J/kg) 
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APPENDIX  C. 


REQUIRED  INPUT  PARAMETERS  FOR  DETERMINING  CAVITY  PENETRATION 


INPUT  PARAMETERS 


1 . A1 , 

2.  B1 , 

3.  C8, 

4.  Rl, 

5.  SI, 

6.  T2 , 


2 

bore  area,  (cm  ) 

projectile  base  position  relative  to  breech  face  of  gun,  (m) 
cavity  initiation  position  relative  to  breech  face  of  gun,  (mj 
cavity  radius,  (m) 
projectile  start  of  motion,  (ms) 

Taylor  coefficient 
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